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Abstract— The analyses presented in this paper are focus on the solution 
of cases I and 3A proposed by the second Aeroelastic Prediction 
Workshop (AePW-2), using an open-source CFD code. The reference 
cases presented by AePW-2 analyze the transonic flow around a 
Benchmark Supercritical Wing (BSCW). AePW-2 Test case 1 consists of a 
forced oscillation problem with Mach number of 0.7 and angle of attack of 
3 deg, while AePW-2 Test case 3A analyzes a flow with Mach number of 
0.85 and angle of attack of 5 deg, being that an unforced and unsteady 
problem. In the study, we simulated both test cases using the software 
SU2, being the results validated by comparison with experimental data 
provided by AePW-2. The results matched with accuracy with the 
experimental data and presented a good response for the analyses of 
AePW-2 test case 3A, proving the software capability of capture the 





physical phenomena involved in this type of flow. 





I. INTRODUCTION 


Computational Fluid Dynamics (CFD) evolved a lot 
during the past two decades. To keep the improving state 
art of CFD, institutions around the world are developing 
workshops, among them, and Aeroelastic Prediction 
Workshop series (AePW) stands out, [1] provides more 
information about AePW. 


The focus of the first edition of the AePW workshop 
series was the solution of unsteady aerodynamics problems 
over three different wing geometry (the Rectangular 
Supercritical Wing, the Benchmark Supercritical Wing 
(BSCW) and High Reynolds Number Aero-Structural 
Dynamics (HIRENASD)). In its second version, AePW 
focused on the analyses of problems involving flutter over 
the BSCW wing. 


Since 2016, all the studies that presented a complete 
solution of AePW-2 test cases used proprietary codes or in- 
house codes, as seen in [2] and [3]. 


More recent studies, like [4], presented the solution of 
the test cases and expanded these, testing the influence of 
parameter variation but these also using in-house codes. 
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However, proprietary and in-house codes present some 
limitations for academia. In this context, open source 
becomes a better option. But nowadays, the full capabilities 
of open-source codes to solve complex flow problems are 
still unrecognized, with just a few papers given an 
overview of this topic. 


Among the possibilities of open-source CFD codes, 
SU2 emerges as a relevant tool for aeroelastic studies since 
it is focused on aeronautics applications, as presented in 
[5]: 

In [6], the developers of SU2 presented more details of 
the software architecture and capabilities to solve the flow 
problem proposed by two different 
configuration test cases. The focus of [6] was to prove the 


full-aircraft 


capability of the software to solve industry-sized problems. 
But for the current study, the principal importance of [6] 
was proving that SU2 was capable of solving transonic 
flow problems over complex geometries since one of the 
test cases validated was the flow over DLR-F6 Transonic 
Aircraft. 


According to [7], the developers of SU2 focused their 
efforts on verifying the capabilities of the software to solve 
computational 


different test cases of interest in 
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aeroelasticity. The study of [7] analyzed flows over NACA 
0012 airfoil, Isogai wing section, BSCW wing, and also 
presenting the benchmark problems solution for fluid- 
structure interaction (FSI). The importance of the research 
of [7] for the current study was the analysis of the BSCW 
wing test cases, which indicates the capabilities of SU2 to 
solve the test cases of AePW-2. 


In a more recent study of SU2 capabilities of solving 
transonic flows, [8] uses SU2 to develop a methodology 
capable of providing the flow response to small-amplitude 
periodic deformations in a structure. This methodology was 
developed using NACA 64A010 airfoil in transonic flow 
conditions and validated by testing it in an Isogai wing 
section and an AGARD 445.6 wing. The results evaluated 
by [8] were accurate when compared with experimental 
data and other numerical simulation results, reinforcing 
SU2 capabilities. 


Verified the SU2 capability of solving transonic flows. 
The current study aims to expand the usage of open-source 
software to solve complex flow problems of interest for 
aeroelastic analysis. The objective proposed was achieved 
by analyzing the SU2’s ability to solve test cases 1 and 3 
presented in AePW-2 and by comparing the results 
obtained numerically with the experimental data provided 
by the workshop. 


II. METHODOLOGY 


AePW-2 uses the Benchmark Supercritical Wing 
(BSCW) for all the analysis proposed, Fig. 1:. presents the 
BSCW geometry view and its cross-section, a SC(2)-0414 
airfoil. This rectangular wing has a chord of 0.4064 m, a 
span of 0.8128 m, a reference area of 0.3303 m?, and a 
moment reference in (0.1219, 0, 0) m. 


The BSCW configuration presents geometric 
simplicity, allowing to set the focus of AePW-2 on flow 


behavior. 


[9] provided the experimental data of wind tunnel 
analysis for test cases | and 3A of AePW-2, being these 
evaluated for a cross-section of the wing, distancing 
0.48768 m from the wing root. Table 1 synthesizes the 
information about the test cases verified in the current 
study, used to test SU2 capabilities. 
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Fig. 1: Benchmark Supercritical Wing (BSCW) geometry 
used by AePW2 (presented in [1]]). 


Table. 1: Test Cases Proposed by AePW-2. 


Case | Case 3A 
Mach Number 0.7 0.85 
(Ma) 
Angle of Attack 3° 5° 
(AoA) 
Fluid R-134a R-134a 
Data type Forced Unforced 
Oscillation Unsteady 
Reynolds Number 4.560- 10° 4.560: 10° 
(Re) 
Freestream 118.0588 m/s 118.0588 m/s 
Velocity (V) 
Speed of Sound (c) 168.6556 m/s 168.6556 m/s 
Temperature (T) 304.2128 K 304.2128 K 
Density (p) 1.1751 kg/m? 1.1751 kg/m? 
Sutherland 243.3722 K 243.3722 K 


Constant (C) 


Reference dynamic 1.1165- 10° Ns/m? 1.1165- 10°° Ns/m? 
viscosity (LUrer ) 


Reference 273 K 273 K 
Temperature (Tref ) 


All the experimental data for the test cases presented in 
Table 1 are from NASA Langley Transonic Dynamics 
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Tunnel (TDT). The test case 3 points to shock-induced 
separated flow in the upper surface and the aft portion of 
the lower surface for Ma = 0.85 and AoA = 5°. 


2.1. MATHEMATICAL MODEL 


Since all the analyzed test cases use R-134a is possible 
to consider the fluid as an ideal gas. Adopting this 
hypothesis is possible to create a correlation between the 
dynamic viscosity (u) and the absolute temperature (T), via 
Sutherland’s law, defined in (1). 


a= ( T ) “Tarp + 5 
= Href] aye 
Tor} THS D 


In all the AePW-2 test cases, the fluid flow is 
considered turbulent. To model the turbulence, we adopted 
the Reynolds-averaged Navier-Stokes equations (RANS). 
With that approach, the governing equations fall on a 
closure problem. To solve this, we used a turbulence 
model. 


Based on the study of [3], was used the Spalart- 
Allmaras Turbulence Model for the analyses of case 1 in 
steady condition and case 3A. The Spalart-Allmaras model 
is a one equation model defined according to [10] by the 
equation (2). 


TẸ 





ar or ss Ea (vy 
At qotu i Bx; 5 = ml — fels — (Corfe - — fa) (= ) 
Me av av av 
= ax; wt ge Ši h ay. Ox; Ox; (2) 
Being the turbulence viscosity (u) defined as: 
Me = PY fri (3) 
Where f yı and y are determined as: 
agg 
X 
fi = FT 
A T Cyl (4) 


2 | 


(5) 
For this turbulence model, the adopted boundary 
conditions are: 


v wall =0 (6) 


0.210438 - væ = Y parisa E 1.294234 + Vo 


Since the Spalart-Allmaras turbulence model is a one 
equation model, it is considerably faster than other models 
with more equations. 


[10] presents the constants and auxiliary relations for 
the Spalart-Allmaras Turbulence Model. 
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For case 1 transient simulation, we considered the 
turbulence model proposed by [11], the shear stress 
transport, or k — œ SST, which is a two equations eddy- 
viscosity model. This formulation consists of a set of 
equations for turbulence kinetic energy and the specific 
dissipation rate equations complemented by the kinematic 
eddy viscosity equation, given by (8), (9), and (10). 








_= pak 
= max(a,w; SF) A 
a( pk. a k, a 
at a Xj ax; 


fe: + soot 








(9) 
(ow) y alow) __ a 
p ðk dw 
i + Ta) + 2(1 — Foa 
Aj w Ox, ax, (10) 


[11] presents more detail about the coefficients and 
auxiliary relations for the k — œ SST turbulence model. 


2. 2. COMPUTATIONAL ANALYSIS 
2.2.1. Mesh 


We generated the mesh using the Ansys Mesh, from 
Ansys academic license, software details can be found in 
[12], and verify the uncertainty due to discretization 
calculating the Grid Convergence Index (GCI), following 
the procedure proposed by [13]. 


For all the meshes developed, we centered the wing 
profile in a semispherical farfield, as can be seen in Fig. 2. 
The figure also presents the boundary conditions adopted 
in the analysis. Table 2 shows the parameters used in the 
mesh generation for cases 1 and 3. 


For case | grid convergence analysis, we developed the 
coarse, intermediary, and fine meshes present respectively: 
156819 elements, 426703 elements, and 1184414 elements. 


The obtained refinement factor was: 1.405 between the 
fine and the intermediary mesh; and 1.396 between the 
intermediary and the coarse mesh. 
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Fig. 2: Mesh developed for test case 1. 


Table. 2: Test Cases Proposed by AePW-2. 


Parameter Case 1 Case 3A 
y+ ] l 
Aspect Ratio 1,2 12 
Number of elements in the boundary 35 35 
layer 
Farfield radius (m) 20 20 


First element height (m) 2,43-10° 2,47-10° 


Following the calculation procedure proposed by [13] 
we estimate uncertainty due to discretization using the GCI 
and obtained the results presented in Table 3. 


Table. 3: Parameters obtained for estimate uncertainty due 
to the discretization of the BSCW wing. 


Refinement factor r21 1,405 

Refinement factor r32 1,396 
Approximate relative error €a21 0,93 % 
Approximate relative error €432 13,12 % 
Extrapolated relative error eex21 0,07 % 
Extrapolated relative error €ex32 1,03% 
Convergence index GCI; 0,085 % 
Convergence index GCI32 1,272 % 


Comparing the parameters presented in Table 3 with 
the exhibit in [13], we saw that the convergence index 
allows the use of the intermediary mesh for all the 
calculations. Based on that result, we developed the meshes 
for case 3A the difference, in this case, was the use of a 
refinement box around the wing, as presented in Fig. 3. 
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With the adoption of a refinement box, we did a local 
refinement in the mesh to capture flow features of pressure 
distribution around the wing. The most dominant feature 
found in the flow is the shock-waves dynamics that should 
occur at the Mach number of 0,85. With this refinement, 
the mesh developed for case 3A had 1768317 elements, 
and the focus of this the upper region of the wing to 
capture the shock-wave dynamics. 


2.2.2. Software 


We used Ansys Mesh from Ansys License of Ansys 
2017 for the mesh generation, [12] presents details about 
this software. 





05 (mm) 0,000 L008 2,000 (rn) 


0,500 1500 


Fig. 3: Mesh developed for test case 3A. 


For the numerical simulation, we used SU2 version 
v6.2.0 Falcon to solve the Navier-Stokes equations. [5] 
presents more detail about the software. 


We evaluated the solution with the following settings: 
Green-Gauss numerical method to compute the gradient; 
FGMRES with ILU preconditioner to solve the linear 
system; JST as flow convective numerical method and 
Scalar Upwind as the turbulent convective numerical 
method. 


For the post-process, we used Paraview 5.7.0. [14] 
provides details about Paraview. 


III. RESULTS AND DISCUSSION 
3.1. Case 1 


In Fig. 4. are presented the results obtained with the 
numerical simulation of test case 1 for the steady flow 
condition. For this test condition, we sampled 76 points 
over the analyzed section and compared them with the 35 
points found in the experimental data provided by [1]. 


As can be seen in Fig. 4., the numerical data almost fit 
with the experimental data provided by AePW-2 for the 
lower surface of the airfoil. 


For the upper surface, numerical and experimental data 
present the same behavior in the Cp curve but diverges in 
magnitude. This divergence in the upper surface occurs 
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because the tetra/prism mesh generated kept some lower 
quality elements in the region. 


Also, Fig. 4. shows that on the trailing edge of the 
wing, the numerical simulation diverges from the 
experimental data. This problem occurs because of the 
sharper edge used in the geometry model. Due to that fact, 
the software couldn’t generate good quality elements, 
leading to an increase in numerical error. 


1,5 
= Numerical Simulation of Case 1A 
—@ Nasa Experimental Data for the Upper Sur- 
face 


~*~ Nasa Experimental Data for the Lower Sur- 
face 





0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1 


x/c 
Fig. 4: Cp plot for numerical and experimental data of test 


case I. 


With the results, we concluded that SU2 could solve the 
steady transonic fluid flows with great accuracy since, in 
Fig. 4., we saw that most of the issues took place due to 
poor quality elements generate in some regions of the 
geometry. 


The major problem found for the analysis was the mesh 
generation. This issue occurs due to SU2 uses meshes in 
SU2, CGNS, and NETCDF_ASCII formats, and just a few 
software develop great quality mesh in these formats. 


During the study, we found that Ansys mesh was the 
only software capable of generating meshes for SU2. We 
also tested Gmsh, but at that time, it didn’t generate proper 
meshes. For this reason, we used Ansys mesh to develop 
all the meshes for the studied test cases. 


For case 1 transient condition, was verified the forced 
oscillation occurring over the BSCW wing. We simulated 
this condition with an oscillation frequency of 10 Hz and 
an angle of 1°. Fig. 5. presents the pressure coefficient 
evaluated with the numerical analysis, and we can compare 
this with the pressure coefficient found by [3] for the same 
test case, exposed in Fig. 6. 


As can be seen in Fig. 5. and Fig 6. the results 
evaluated by the authors keep the same behavior as the 
results evaluated by [3]. 


The magnitude of the peak curvature is analogous to the 
one found in [3]. However, the curvature found by [3] 
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presents two peaks, while the curves obtained by the 
authors present a single peak. Again the pressure 
coefficient next to the trailing edge was poorly represented 
in comparison with the found by [3]. 





0 0.2 0.4 0.6 0.8 1 
x/C 


Fig. 5: Cp coefficients obtained by the authors for test case 
I transient condition. 





“9 0.2 0.4 0.6 0.8 1 


Fig. 6: Cp coefficients obtained by [3] for test case 1 
transient condition. 


Another way to see the behavior of SU2 is to plot the 
results in the frequency spectrum. AePW-2 presents the 
frequency response at 10 Hz for the sensors applied in the 
experimental tests. We can see a comparison between this 
response and the computational responses obtained by SU2 
in Fig. 7. and Fig 8. 


In Fig. 7. and Fig 8., we can see that the values 
obtained by SU2 are similar to the experimental evaluated 
by AePW-2, keeping the same shape and same peaks at 
upper and lower surfaces. 


3.2. Case 3A 


Since case 3A consists of an unsteady problem, it was 
necessary to adopt a time step for developing the 
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interactions over time. For the analysis, we used a time step 
of At = 10* s. Fig. 9. presents the results obtained for the 
SA model and Fig. 10. for the k-—w SST model. 


—=FRF-Magnetude Su2 {Lower} 


——FRE-Magnetude AePW2(Lower) 


0,06 


0,04 


0,02 





0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1 


Fig. 7: Comparison between the magnitude frequency 
response at 10 Hz for the lower surface. 


—=FRFE-Magnetude Su2 (Upper) 
——FRE-Magnetude AePW2(Upper) 


0,25 


0,1 


0,05 





Fig. 8: Comparison between the magnitude frequency 
response at 10 Hz for the upper surface. 


As presented in Fig. 9. and Fig. 10., the numerical 
results almost fit with the experimental data for this case. 
The difference found stays on the transition of the Cp that 
occurs next to x/c = 0.16, where the experiments present an 
abrupt fall of the Cp, while the numerical results exhibit a 
smooth transition. 


Comparing case 3A and case 1 results, it is possible to 
see that the first presented more accuracy due to the mesh 
used. 


Since case | consists of a flow with a low Reynolds 
number, and the problem occurs at a steady-state, the mesh 
for this case was coarser than case 3A mesh due to it 
doesn’t use the refinement box. These simplifications into 
the mesh reduce the computational cost but sacrifice part of 
the solution’s accuracy. 


For case 3A, since the problem involves capture the 
shock wave dynamics over the wing was necessary to 
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adopt local refinement techniques in the mesh generation. 
Due to the local refinement, we minimized the trailing edge 
problem found in case 1 and got a more accurate solution. 





— Numerical Simulation of Case 3A 


—~@ Nasa Experimental Data for the Upper | 
Surface 


1 =- Nasa Experimental Data for the Lower 
Surface 





0 0,1 0,2 0,3 0,4 


x/c 
Fig 9: Comparison between Cp plot for numerical and 


experimental data of test case 3A using SA model. 


= Numerical Simulation of Case 1A 
=—@ Nasa Experimental Data for the Upper Sur- 
face 


~- Nasa Experimental Data for the Lower Sur- 
face 





0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9 1 
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Fig 10: Comparison between Cp plot for numerical and 
experimental data of test case 3A using k—q@ SST. 


Another detail noticed is the difference evaluated by the 
turbulence models. While the SA model captured the Cp 
variation over time, as seen in Fig. 9., the k — œ SST 
wasn’t capable of that, as presented in Fig. 10. 


Also, Fig. 9. and Fig. 10. presents that despite both 
turbulence models represent the behavior of the flow over 
the wing adequately, but none captured the discontinuity 
presented by the shock wave. 


IV. CONCLUSION 


After all the analyses, we confirmed the capability of 
SU2 to solve transonic problems. 


During the study, the principal limitation found was the 
generation of a proper mesh. Since SU2 native format is 
.su2, our first attempt was to use open-source mesh 
generators capable of generating meshes in this format. 
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None of the .su2 Open-source mesh generators tested 
generated meshes that provided good results for SU2. 


Due to that, during the study were necessary to use 
another mesh format. In this case, was used the CGNS 
format, being the meshes generate by Ansys Mesh. 


The results also present that the generated mesh 
impacts the accuracy of the simulation. Since a more 
refined mesh, like the one used for the numerical 
simulation of case 3A, was more accurate when compared 
with the coarse mesh generated for case 1, even 
considering that complexity of case 3A greater than case 1. 
This result also shows the importance of local refinement 
for unstructured meshes. 


The analysis of case 3A presents that SU2 was capable 
of capture the shock wave dynamics. Also, the numerical 
results almost fit with the experimental data provided by 
the workshop AePW-2. 


As observed in Fig. 9. and Fig. 10., the major problem 
found for the analysis was the capture of the abruptly falls 
off the Cp over the upper surface of the BSCW wing since 
the numerical simulation presents a smooth transition 
between the Cp curve while the experimental data shows a 
more abruptly fall. 
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